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We present here the electronic structure and optical properties of InGaAs quantum wells with barrier doped 
with Manganese. We calculated the electronic states and optical emission within the envelope function and 
effective mass approximations using the spin-density functional theory in the presence of an external magnetic 
field. We observe magneto-oscillations of the Landau levels at low-magnetic fields (J? < 5 T) that are dominated 
by the magnetic interaction between holes spin and Mn spin, while at high magnetic fields the spin-polarization 
of the hole gas is the dominant effect. Our results also show that a gate voltage alter significantly the magneto¬ 
oscillations of the emission energy and may be an external control parameter for the magnetic properties of the 
system. Finally, we discuss the influence of the Landau Levels oscillations in the emission spectra and compare 
with available experimental. 


PACS numbers: 73.21.Fg, 78.67.De, 78.55.Cr, 75.50.Pp 

I. INTRODUCTION 

The research on magnetic semiconductors has attracted 
much attention for more than two decades'"^. The most 
investigated material is the (Ga,Mn)As system, where it 
has been observed ferromagnetic phase with Curie temper¬ 
atures (Tc) reaching 190 K for samples with Mn concen¬ 
tration of ^ 10%.^’^ Others (III,Mn)V materials, such as 
(In,Mn)As, (Ga,Mn)Sb, (In,Mn)Sb, have also shown ferro¬ 
magnetic phase®. In these materials, Mn acts as both an ac¬ 
ceptor and a magnetic impurity and its ferromagnetism is me¬ 
diated by the interaction between holes and Mn spins^’’. As 
the magnetic interactions are mediated by charged carriers, the 
control of the magnetic properties can be achieved by electri¬ 
cal and optical means*'^®. 

Much of the research efforts have been concentrated on 
(III,Mn)V bulk or in its heterostructures where both Mn and 
holes are in the same spatial region. These structures allow 
a high ferromagnetic Curie temperature due to the strong in¬ 
teraction between the hole gas and the Mn ions. However, 
the hole gas is strongly scattered by the Mn ions, reducing its 
mobility and the optical quality. To overcome this difficulty, 
the (Ga,Mn)As layers were grown in the presence of a quan¬ 
tum well, as for example in a GaAs-(In,Ga)As-(Ga,Mn)As se¬ 
quence. In this situation, the hole gas is located in the quantum 
well, separated from the Mn ions which need to be controlled 
in order to maintain a certain level of overlap between the 
holes and the Mn ions to assure the magnetic properties. Re¬ 
cently, such heterostructures have been investigated by means 
of transport"’*^ and optical experimentsThe results sug¬ 
gest that the interaction between holes and Mn is determinant 
for these systems properties. 

Gazoto et al. investigated (In,Ga)As QWs with GaAs bar¬ 
riers (5-doped with both carbon and Mn in alternate sides of the 
QW. The samples were 5-doped with Mn in order to increase 
the Mn doping concentration beyond the solubility limit. The 
presence of the 5-doped C layer in the other side of the QW 
aimed to increase the hole gas concentration and, with that, 
to increase the magnetic effects. They observed that the cir¬ 
cularly polarized magneto-photoluminescence presents strong 


oscillations with the magnetic field and they are more pro¬ 
nounced in the samples with higher Mn concentration. These 
oscillations were attributed to the Landau level filling fac¬ 
tor. Magneto-oscillations of the circularly polarized emission 
were also observed in both two-dimensional electron (2DEG) 
and hole gas (2DHG)'*“^^®. The origin of these oscillations is 
in the many-body effects of both the two-dimensional gas^'’^^ 
and in the optical recombination process^^“^®. The oscilla¬ 
tions of the transition energies observed in Ref. 14 are much 
stronger than those observed in other high quality 2DHGs^®’^®. 
These oscillations were correlated to the presence of Mn spins 
in the heterostructure. Several other experiments such as low- 
magnetic field circularly polarized, photoluminescence exci¬ 
tation and time-resolved photoluminescence*^’^® and trans¬ 
port measurements^^, were performed in similar heterostruc¬ 
tures. They all showed significant magnetic interaction corre¬ 
lated to the presence of the hole gas. 

Here, we present the results of a calculation for the elec¬ 
tronic structure and the emission energy of (In,Ga)As QWs 
with GaAs barriers 5-doped with Mn and C. The electronic 
states of the heterostructure are calculated using the spin- 
density function theory (SDFT)^®^^^ within the envelope func¬ 
tion and effective mass approximations^^. The Mn-hole 
spin interaction is described by the Zener kinetic exchange 
theory^’^’^^. We aim to provide useful information to under¬ 
stand this system and to obtain a microscopic interpretation of 
the observed effects. We compare our results with the experi¬ 
mental results from Gazoto et al. and suggest an interpreta¬ 
tion for the observed magneto-oscillations they observe in the 
emission spectra. 

The paper is organized as follows. In Section II we present 
the model used to obtain the electronic states of the structure 
and the emission energies. In Sec. Ill we discuss in details the 
theoretical results focusing on the samples used by Gazoto et 
Finally, in Sec. IV we present our concluding remarks. 


II. MODEL 

The model heterostructure investigated is illustrated in 
Fig. 1(a) and it is based on the system studied by Gazoto et 
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It is composed by a 500 nm GaAs buffer layer, a car¬ 
bon (C) 5-doping layer, followed by 10 nm GaAs spacer, a 10 
nm Ing lyGag 83 AS layer, a GaAs spacer Lg, a Mn 5-doping 
layer with a concentration xmu given in percentage of mono- 
layers (MLs), and finally a 60 nm GaAs top layer. The dis¬ 
tance Lg between the QW and the Mn doping layer and the 
Mn concentration will be considered as free parameters in our 
model to study the effects of the Mn interaction. All the others 
parameters are hxed. As a consequence of the Mn and C 5- 
doping and the thermodynamic equilibrium, the QW presents 
a 2DHG. We consider the system under an external magnetic 
held applied along the growth direction, here the z-direction. 


body problem in a set of non-interacting equations, which are 
solved self-consistently. 

We approximate the sample as being homogeneous in the 
plane. The z-direction and the in-plane (x, y) directions are 
therefore not coupled. The Hamiltonian can be written as 

jjhh{lh) _ jjhh{lh) _j_ ^hh{lh) 

where the hrst term is the z-part of the Hamiltonian and the 
second term is the in-plane Hamiltonian. The xy-part of the 
Hamiltonian is responsible for the formation of the Landau 
Levels (LLs)^®. The Hamiltonian in the z-direction is 
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FIG. 1. (Color online) Schematic representation of the investigated 
heterostrucutre and its self-consistent potential profile and wave- 
functions. 


We used the envelope function and effective mass approxi¬ 
mations to describe the electronic states of the QW^^. The al¬ 
loy is treated in the Virtual Crystal Approximation. In general 
to obtain the valence band (VB) states it is necessary to use 
the six-band Luttinger-Kohn Hamiltonian^^ that describes the 
heavy-hole (hh), light-hole (Ih) and split-off {so) bands and 
the coupling among them as well as the spin-orbit effects^. 
However, in our system, the Ino iyGaossAs layer is under 
compressed stress which splits the hh and the Ih bands by a 
value of the order of 50 meV. The Ih band is actually a type-II 
or marginally type-I heterostructure, depending on the param¬ 
eters chosen to describe the structure. This results in a near 
parabolic dispersion for the hh QW ground state for the ener¬ 
gies of interest to study the optical emission. We therefore 
simplify our approximation considering a simple parabolic 
dispersion and tested it against a full Luttinger Hamiltonian 
calculation at zero magnetic field. The comparison is very 
good for the ground state responsible for the optical emission. 
As we go for the excited states, some discrepancies may be 
present but they do not alter the main results we show here. 

The many-body effects of the 2DHG are considered within 
the SDFT^°’^^. This allows us to calculate the ground-state 
properties of the 2DHG including exchange and correlation 
(XC) effects in the presence of a spin-dependent potential. 
We employed the Kohn-Sham minimization scheme^® to ob¬ 
tain the electronic structure. This procedure maps the many- 
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The first term is the kinetic energy, the second is the het¬ 
erostructure potential, the third and fourth terms are the 
Hartree and XC potentials, respectively. The fifth is the Zee- 
man contribution, and lastly, the hole-Mn (h-Mn) coupling. 
We will now describe each of the Hamiltonian terms in de¬ 
tails. 

Heterostructure potential - is the structural po¬ 

tential which is built up from the band gap difference and band 
alignment between Gao.salno.irAs and GaAs layers plus the 
strain effects. The energy gaps, at low-temperature, are 1.264 
eV^’ and 1.519 eV^^ for Gao.ssIng.irAs and GaAs, respec¬ 
tively. The band alignment is type-I. As the QW is strained, 
its total band offset contains both the band gap alignment and 
strain contributions. We first define a band offset without 
strain^^. In this case, we assume a VB (CB) offset. Ay (Ac) 
of 15% (85%) of the energy gap difference. 

Our system is dominated by the GaAs layers and we as¬ 
sume that the whole structure presents the GaAs lattice pa¬ 
rameter. This generates a compressive biaxial strain in the 
Gao. 83 Ino. 17 As layer, which alters the QW band offset^®’"^®. If 
we neglect the so band, T 7 , the effect of the compressive bi¬ 
axial strain is manifested in an hydrostatic term {5Eh), which 
increases the gap, and a shear deformation {6Eg), which splits 
the hh and Ih bands'^®. The deformation potential for the hy¬ 
drostatic term can be split in a contribution to the conduction 
band, SE^^ and one to the valence band, SE^^. The final 
result is a QW and a barrier potential for the hh and Ih bands 
in the Gao. 83 Ino. 17 As, respectively. The hh and Ih potential 
are given by 


-rrhh j Ih 
^QW 


Av - SE^^ ± SEg, 


(3) 


where 


5L//^ “F tlyy “F Cjjjj;) 

SEg = by(^€xx “F t:yy ^CzZ^)- 

For compressive biaxial strain, the strain components are 

given by Cxx — — (^GaA.s tZCalnAs)/tiliiGaAs ^ 0? 

and = - 2 C' 2 iea:a:/C'ii, where ocaAs and ocainAs are 
the GaAs and Gao. 83 Ino. 17 As lattice parameters, respectively. 
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avB(CB) and by are the deformation potentials. Cu and C 21 
are elastic stiffness constants^®“‘^°. 

We can now turn our attention to the insertion of the Mn 
in the GaAs. The Gai_a;j^^Mna;j^^^As layer has a larger lat¬ 
tice parameter than GaAs and follows a similar analysis as 
for the Ino.irGao.saAs layer regarding the strain effects. One 
additional difficulty is that the Mn strongly diffuses towards 
the surface in the GaAs as it was shown from secondary ion 
mass spectroscopy (SIMS)"^'^^. Based on the SIMS results, 
instead of considering a Mn (5-doping layer, we assumed that 
the total Mn is distributed over many GaAs layers, forming a 
Gai_a;„_jMna;jy^As alloy. To take this effect into account, we 
assume that the system is homogeneous in the x-y directions 
and only variations of the Mn concentration in the z-direction 
are considered. We construct a distribution function that takes 
in account the Mn diffusion in the GaAs^*^^. The SIMS 
results show that the Mn diffuses following an approximate 
gaussian-like function in both directions of the heterostruc¬ 
ture. However, it diffuses more strongly in the direction of the 
surface than towards the QW. We consider therefore a double- 
gaussian distribution function, as defined below 

f{z - L,) = foexp{-[{z - Ls)/A(z)]2}, (4) 


where Lg is the gaussian center, that is, the nominal (5-Mn 
doping position. A(z) is the average width of the gaussian 
that describes the Mn diffusion 


A(z) 


d \f z < Lg, 
D if z > Lg, 


(5) 


and /o is the normalization constant. We considered D = 2 
nm and d = 1 nm, which are compatible with the SIMS 
results^'^^. The Gai_ 2 :^j„Mna;„_jAs layers also shows a dif¬ 
ferent gap than GaAs, and an intrinsic band offset should be 
present. The main contribution for the Gai_a;j^j_^Mna;^j_^As 
band offset, however, has origin in the sp-d interaction, which 
is discussed below. We therefore will neglect the intrinsic 
band offset. 

Finally, we consider the S-C layer. C act as an acceptor and 
its growth is well controlled and does not significantly diffuse. 
We assume an homogeneously distributed 5 A doped region 
in a continuous approximation. The heterostructure potential 
finally can be written as 


- ^Qw/4) 

- T <5i?«T)aGaMnAs/(^ - Lg) 


with Q{x) being the Heaviside function, and Lqw the QW 
width, that is, the Ino. 17 Gao. 83 As layer. 

Hartree potential - The third term of is the Hartree po¬ 
tential, which is obtained by solving Poisson’s equation 

- NcQ{z - zl)e{z° - z) ( 6 ) 

-PMnf[z - Lg)] , 


where e is the GaAs dielectric constant. p{z) is the total hole 
density, which at zero temperature (T = 0 K) is given by 

if S = OT, 

b,i 

P(^) = S Jl to'’ ^ . 

E ^(") 1^ /-E 4". if i? A OT. 

b,i 

< Ll-,Tz 

(7) 

where 



eB 1 
Stt/i v^r 



2 r2 


( 8 ) 


The first (second) line describes the hole density in absence 
(presence) of an external magnetic field, ^^{z) are the VB 
envelope functions, and Aj^ = Ep — , where Ep is 

the Fermi level, and 7?^ the eigenvalues of H^{b = hh, Ih). 

= E’P^^ + E^ is the total subband energy of the &-hole 

with spin Tz in the n-th LL of the Ath subband. F = Fo vA? is 
the LL broadening related to the 2DHG mobility'*^'^^ in our 
calculations, we considered Fq as a parameter. 

In the absence of a gate voltage, we assume that the Fermi 
level, Ep, is pinned at the surface states, that is, in the middle 
of the gap at the surface, Ep = EGaAs 126 AS _ ^i^q 

consider the case of an applied gate voltage, Vg, which allows 
us to change the Fermi level position. In this case we have 
Ep = Ef^^^j2 - Vg. 

Exchange-correlation potential - The forth term of Eq. (2) 
is the XC potential. Here we use the Vosko, Wilk, Nusair 

(VWN) parametrization^®47 fQj- ji^g local-spin-density approx¬ 
imation. vxc{z) depends on both the hole density p{z) = 
p^{z) -\-pi{z) and the hole gas magnetization ^(z) = p^(z) — 
Pl{z). Again, we consider that the density in the plane is ho¬ 
mogeneous, and hence, the hole density depends only on the 
z-coordinate. 

Zeeman potential - The fifth term in Hz is the Zeeman in¬ 
teraction between hole spins and the external magnetic field. 
g* is the hole effective g-factor and p ,b the Bohr magneton. 
The lih and Ih spins are = ±3/2 and r]'’ = ±1/2, re¬ 
spectively. 

p-d potential - The last term of the Hamiltonian is the p- 
d interaction between holes and Mn spins. This term has its 
origin in the interaction between VB states with the d-orbitals 
of the Mn impurity'**^’"'^. The presence of the hole gas is de¬ 
scribed via Zener kinetic-exchange modeA’^’^^. 

We assumed the magnetization vector aligned along the 
z-direction. This is a valid approximation since even a 
small magnetic field aligns the magnetization in its direction^" 
against the easy-axis'^^’^'*. 

The final expression is written as 

VpT'^\z) = -^Nol3xyffMT^^^^^'>BM{y)aMnf{z-Lg), 

(9) 

where Ao/3 is the p-d exchange constant of the spin interaction 
between Mn’s and holes. Xeff is the effective concentration 
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of Mn spins (see below). M = 5/2 is the Mn spin. BMiu) is 
the Brillouin function.^' Its argument is given by 


y = 


yMny-sMB 

ksT 


JpdM 

2 kBT 


Ls)dz, 


( 10 ) 


the first term is due to the interaction of Mn spin with the ex¬ 
ternal magnetic field, where gMn is the Mn g-factor, ks is the 
Boltzmann constant, and T is the temperature. The second 
term is the antiferromagnetic interaction between holes spins 
with Mn spins, which is responsible for the ferromagnetic in¬ 
teraction of Mn’s spins^’^’^^. Jpd = P/Nq is the p-d exchange 
constant. Nq is the cation concentration. 

This antiferromagnetic interaction depends on the 2DHG 
magnetization and the overlap between hole and Mn ions. 
If Mn ions and holes are equally homogeneously distribute, 
the distribution function aMnf{z — Lg) is replaced by a uni¬ 
tary constant and we recover the well-known results from 
(Ga,Mn)As bulk^. Here, they are non-homogeneous and the 
second term in the function y depends strongly on the struc¬ 
tural parameters, namely, the Mn position in the GaAs and the 
holes states. 

Equations 2, 6, and 7 are solved self-consistently. The 
Schrodinger equation is solved via split-operator method^^. 

Conduction band states - The CB states are calculated 
within the same approximations. The in-plane Hamiltonian 
gives the LLs for the CB. The electron Hamiltonian in the z- 
direction is 


correlation energy. We will neglect this term in our cal¬ 
culations. Within our description the electron-hole correla¬ 
tion potential should not present a dependence with the spin- 
polarization of the hole gas, but only on the total hole gas 
density. Therefore, the main effect of this contribution will 
depend on the total charge transfer between the hole reser¬ 
voirs of the heterostructure which plays a minor effect in most 
of our results (see next Section). 

Optical transitions - The transition energy is calculated as 
the energy difference between the electron and hole eigen¬ 
states. Our focus is in the circularly polarized emission. The 
right circularly polarized ((t+) light is given by the recombina¬ 
tion of a spin-down electron with a spin-up hh, while the left 
circularly polarized (ct_) light is the recombination of a spin- 
up electron with a spin-down hh. The recombination energies 
are given by 


^fT+(cr_) _ jpe 

^Tot - 


El 


■E: 


hh 

44) 


■E. 


hh 


The emission is allowed only if the electron and hole states 
are in the same LL (n = m). To enhance the magnetic field 
effects on the transition energies, we subtract the transition 
energy at zero magnetic field from E'^'^^'^~\b), redefining 
the transition energy as 


= E^^J'^-\b) - e^^J'^-Ho). (14) 


We define the non-linear energy shift as 


d? 

HI = - - —^ + vl^^{z) - VB.{z) + gePBCTzB ( 11 ) 

^iiLq CLZ 

+ + t'c(z) 

where irie is the electron effective mass in the CB. The second 
term is the CB heterostructure potential, which can be written 
as 

= (Ac - SE^^)Q{z - Lqw/2)0(Lqw/ 2 - z) 

“'^^°Mn“GaMnAs/(^-is)- (12) 

The strain contribution in the CB offset is limited to the hy¬ 
drostatic term (6Ef^^). The third term is the Hartree potential 
defined in Eq. 6. The fourth term is the Zeeman interaction, 
where pe is the electron g-factor and = ±1/2 is the elec¬ 
tron spin. In the fifth term we have the s-d interaction between 
the electron’s and Mn’s spins, which is written as 

^sdi^) = NoaXeffMa^BM{y)aMnf{z - ±s), (13) 

where A/) a is the s-d exchange constant between electron 
and Mn spin. The last term is the effect of the correla¬ 
tion potential on the conduction band due to the presence 
of the hole gas. Since electrons and holes are treated as 
different particles, there is no exchange contribution for the 
CB. This potential was parametrized for the case of a spin- 
unpolarized 2DHG'^^’^^ and gives an important contribution 
for the band-gap renormalization observed in the optical spec¬ 
trum of modulated-doped QWs^^. Erom our knowledge there 
is no parametrization for the spin-dependent electron-hole 


AE^+^'"-\b) = E^+‘''"-\b) - E1H<^-\b) - El{B) 

(15) 

where we subtracted all the linear terms in B in order to mag¬ 
nify the non-lienar effects in the transition energy. 

Parameters - The holes in the heterostructure are provided 
by both the C and the Mn doping. We will consider the effec¬ 
tive Mn concentrations parameters that better agree with those 
from [14]. The C concentration A), is fixed for all systems we 
investigate here. We consider A), = 13.35 x 10^® cm“® which 
is the value obtained by fitting the measured hole concentra¬ 
tion in the QW for a sample without Mn and comparing with 
our calculations'^. The nominal concentration of Mn, xmu, is 
known from the growth process. However, it does not provide 
the real hole density, since Mn can be either a substitutional 
impurity or an interstitial one^. In the first case Mn replaces 
Ga, and provides one hole to the system, while at the inter¬ 
stitial position it is a double donor, and gives two electrons. 
Therefore, there is a self-compensation of holes by the elec¬ 
trons, and the total density of holes provided by the Mn is 
given by pMn = xs — 2xi, where xs and xi are the concen¬ 
tration of substitutional and interstitial Mn, respectively^. Eur- 
thermore, because of the attractive Coulomb interaction, the 
interstitial Mn ions tends to be near to the substitutional ones 
presenting an antiferromagnetic coupling, which reduces the 
net Mn spins^. The effective Mn spin concentration is given 
by Xef f = Xs — xi^'^ . In our model, we only describe the un¬ 
compensated substitutional Mn, pMn, with the effective spin 
concentration, Xe//. We do not have direct access to xs and 
X/. These values are strongly dependent on sample growing 
conditions. We extract pMn from our calculations by fitting 
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the theoretical value pqw with its experimental value, which 
was estimated from Shubnikov-de-Hass and Stoke shift mea¬ 
surements. This allows to determine PMn for each sample 
which is used as fixed parameters for the remaining calcula¬ 
tions. The other parameters are described in Table I. 

TABLE I. Parameters used in the self-consistent calculation. The 
Ino. 17 Gao. 83 As strain parameters are linear interpolation between 
GaAs and InAs parameters. The parameters were extracted from the 
refs. 7, 37-39, 55-57. 


Parameters 

GaAs 

InAs 

71 

6.98 

20.0 

72^** 

2.06 

8.5 

Cii (10^° Pal^** 

12.21 

8.329 

C 21 (10^° Pa)^** 

5.66 

4.526 

aL 

5.65325 

6.0583 

a^b (eVl^** 

-7.17 

-5.08 

acb (eV)^* 

-1.16 

-1.0 

K (eV)^* 

-2.0 

-1.8 


0.11 



-2.9 



-2.3 


Egixin) (eV)” 

1.519 - L583a:/„ - 1 - 475a;f„ 


Avb (eV)^* 

0A5Eg{xin) 


Noa (eV) 

0.2 


Nofi (eV)''“ 

1.2 


Jpd (meV nm^)^ 

54 


^Mn 

0.4 


^eff 

0.13 


PMn (10^^ cm“^) 

9.92 


Pqw (10“ cm“Q 

5.2 



III. RESULTS 

We consider a sample with XMn = 0.4 ML and Lg = 1 
nm as our case study. The LLs broadening (Tq) is fixed for all 
investigated heterostructures independently of xmu- To is re¬ 
lated to the 2DHG mobility, which is approximately ~ 2 x 10^ 
cm^Ws at 77 K in the QW for the samples investigated in 
Ref 14. This implies in Tq « 1.8 meV The ex¬ 

perimental photoluminescence is performed at 2 K. At lower 
temperatures the mobility increases and this leads to smaller 
values for the LLs broadening. In our calculation we consid¬ 
ered Lq = 0.25 meV Fig. 1(b) illustrates the self- 

consistent profile potential for the hh and Ih VB. The wave- 
functions in the VB represent the occupied /z/i-subbands. We 
note that there is a 2DHG in the carbon doped layer. In the Mn 
layer there is no occupied state. Therefore, the only hole gas 
interacting with the Mn ions has its origin in the QW states. 
We assume that the Fermi level is pinned at the surface states 
which creates a high electric field between the surface and the 
QW as it can be observed in Fig. 1(b). As a consequence of 


this strong potential anisotropy and the different doping at the 
barriers, the QW CB ground state (not shown here) is pushed 
towards the Mn doped barrier. These features are common for 
most of the case discussed here. 

First we examine the evolution of the CB LLs as a function 
of the magnetic field. In Fig. 2(a) we plot the LLs associated 
to the CB QW ground state, g.- At low-magnetic fields, 
more precisely at 1.5 T, there is a crossing of the n=0 CB 
LLs. This crossing originates from the interplay between the 
Zeeman energy and the s-d interaction. It is, however, a too 
small effect to be observed experimentally. At high-magnetic 
fields, B > 5 T, we observe oscillations in the CB LLs, which 
are associated to the crossing of the QW hole LLs with the 
Fermi energy. These oscillations did not change the CB LLs 
dependence with the magnetic field and they show the same 
behavior for both spins. This is expected because the electron 
spin does not interact with the hole spin via exchange inter¬ 
action, and also the s-d exchange coupling cause only a rigid 
energy shift at large magnetic field. Therefore, the observed 
oscillations have their origin in the charge transfer between the 
C layer and the QW, as it will be discussed bellow. We remind 
that the electron-hole correlation was not taken into account in 
our model. However, we do not expect a qualitative influence 
from this term on the spin-properties of the QW hole gas. It 
would certainly lower the energy levels and enhance the os¬ 
cillations since correlation potential should be dependent on 
the total hole density. Its effect is associated to the oscilla¬ 
tions in the total hole concentration. As we observed, this is a 
minor effect [see Fig. 3(a)]. The energy transitions calculated 
within this approximation should give a reliable behavior with 
the magnetic held. 

Next we examine the VB electronic structure. Fig. 2(b) 
shows the LLs fan diagram. In order of increasing energy 
we have at zero magnetic held and hh'^ subbands, 

where the upper index, n and refers to the subband order, 
LL and /z/z-spin, respectively. The first and second subbands 
are located in the QW and in the carbon layer, respectively, 
as depicted in Fig. 1(b). There is also a marginally occupied 
light-hole subband in the carbon layer, which is not shown in 
Fig. 1(b) and does not play any significant role in the subse¬ 
quent analysis. We observe a low- and high-magnetic held 
regimes for the Landau Levels magnetic-field dependence. 

We first consider the low-magnetic field regime, i.e., B < 5 
T. Let us concentrate on hh^ ^ and hh^ ^ states. They are the 
more important states for the emission spectra. In particular, 
for B < 2 T, hh\ ^ state shows an almost flat dependence 
with the magnetic held while hh^ ^ state shows a stronger 
one. The origin of this behavior is three fold; the Zeeman 
effect, p-d interaction, and the spin-polarization of the hole 
gas, which is enhanced by the hole-hole exchange energy. For 
low magnetic fields the spin-up and spin-down hole densities 
are nearly the same, and hence, the 2DHG spin-polarization 
does not influence its magnetic field dependence (see discus¬ 
sion below for high-magnetic fields, where this contribution is 
important). Also, in this regime, the Zeeman energy is small 
and its effect is to split the /z/z-spin states by the same amount 
of energy. Therefore, the hole-Mn spin interaction is the dom¬ 
inant effect at low magnetic fields. As the magnetic held in- 
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creases, the Mn spins are aligned and the Brillouin function 
saturates. The net result is a sizable spin-splitting originated 
by the p-d coupling that remains at higher magnetic fields. 
The flat behavior of hhQ ^ originates from the competition be¬ 
tween p-d interaction and diamagnetic shift of the LLs, which 
tends to cancel out, while for the hhg ^ they add up. This was 
verified by turning off p-d interaction (not shown here) in our 
model. 
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FIG. 2. (Color online) (a) Electron Landau levels fan diagram. The 
solid and dash-dot lines represents the .j. and ^ LLs, respec- 
tively.(b) Fleavy-hole Landau Level as a function of the magnetic 
field. Solid, dash-dot, dash-three-dots, dash two-dots, and dashline 
represents ™d Ef, respectively, (c) Tran¬ 

sition energy (left i/-axis) and non-linear energy shift (right y-axis) 
as a function of the magnetic field. Solid and dashdot line are the 
transition energies for and (t_ polarization. The non-linear ener¬ 
gies AEa^ and AEcr_ are represented hy open and closed circles, 
respectively. 


We turn our attention now to the high-magnetic field 
regime, i.e., B > 5 T [see Fig. 2(b)]. We observe that all 
LLs and the Fermi level oscillate as a function of the mag¬ 
netic field. These oscillations have origin in both the crossing 
of the LLs with Ep, and the charge transfer between the hole 
gas in the QW and the carbon layer. However, we observe that 
hh^ ^ and hh^ ^ oscillate in the same or opposite directions, 
depending on which spin hole level is crossing Ep. On the 
other way, the /i/ig and ^ levels roughly following the 
Ep oscillations. Fig. 3(a) shows the QW and the C 2DHG 
concentration as a function of the magnetic field. It gives us 
a measure of the charge transfer between the QW and the C 
layer as a function of the LL filling factor. Fig. 3(b) shows the 
LLs occupation as a function of the magnetic field. Most of 
the LLs show a similar behavior. Their hole concentration in¬ 
creases linearly with the magnetic field, as it is expected from 
the LL degeneracy. As one LL crosses the Fermi level, it starts 
to be depopulated. If there was no broadening, this should be 
an abrupt decrease. In our case, the broadening makes the de¬ 
population of the LL to last a finite range of magnetic field, 
but with a nearly linear decrease. This behavior is consistent 


for all LLs associated to QW subbands. The C layer LLs show 
a different behavior. As they start to be depopulated, they do 
not follow a linear behavior. Actually, this behavior is associ¬ 
ated with the C layer LLs states that roughly follow in energy 
the Fermi level, as if they were partially pinned on it. 
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FIG. 3. (Color online) (a) Total (pr), quantum well (pqw), and car¬ 
bon (pc) two-dimensional hole concentration, and (b) Landau level 
two-dimensional density as a function of the magnetic field. 


Taking that in consideration let us look now in more de¬ 
tail in the QW LLs oscillations and their correlation with the 
LL filling factor. Here again we focus on the analyses of the 
lowest energy QW LLs, i.e., hh^ and hh^ ^ since they are 
responsible for the optical emission. At B ^ 5 T the states 
at the QW that are fully occupied are , i = 0,1 and 
Tz =t! We have in this situation a spin-unpolarized hole 
gas in the QW. Other LLs from states at the C layer are also 
occupied but they do not affect the results we discuss here. At 
B ~ 5 T, as depicted in Fig. 3(b), hhl ^ starts to be depopu¬ 
lated, spin-polarizing the QW hole gas. It depopulates entirely 
at i? ^ 7 T when the hole gas in the QW becomes partially 
spin-polarized. In this interval, 5-7 T, /i/ij and hh^ ^ have a 
significant difference in their magnetic field dependence, with 
hh^ ^ energy increasing strongly with the magnetic field while 
hhQ shows a weak dependence with it. At i? ~ 7 T hhl t 
starts to be depopulated and the QW hole gas starts to de¬ 
crease its spin-polarization until 5 ^ 11 T when hhl^^ is 
completely emptied and the QW hole gas is spin-unpolarized. 
In this interval of magnetic field, both L/iJ and hhQ ^ states 
show similar weak magnetic field dependence. At H ^ 11 T 
hhQ ^ starts to be depopulate and the hole gas is again spin- 
polarized. The magnetic field dependence of L/iJ and /i/ij ^ 
states again differ significantly, repeating the previous pat¬ 
tern. The spin-polarization and spin-unpolarization manifests 
itself by an oscillation in the LL dependence with the mag¬ 
netic field. These oscillations take place each time a QW LL 
is emptied and a new one starts to be depopulated. 

It should be observed that during this range of magnetic 
fields 5-11 T, hh and Ih LLs associated to the C layer are also 
changing their occupation in relation to their maximum occu¬ 
pation but that does not affect the L/iJ and hhQ ^ magnetic 
field dependence. Actually, this picture is confirmed by the 
charge transfer between the C layer and the QW. Fig. 3(a) 
shows that the charge transfer oscillates following the QW LL 
filling factor. This charge transfer has its origin in the thermo- 
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dynamic equilibrium. However, it is not significant due to the 
weak overlap between hh^ and hh"^ wave-functions. Clearly, 
the most important effect is the spin-polarization of the hole 
gas in the QW, a consequence of the charge transfer between 
LLs within the QW. 

We now turn our attention to the consequences of these ef¬ 
fects on the optical emission. Fig. 2(c) illustrates the fun¬ 
damental energy transition shift, and the non¬ 
linear energy shift as a function of the magnetic 

field. We focus our attention on the non-linear energy shifts 
where the magnetic effects are more clearly displayed. We 
first observe at low magnetic fields a strong non-linear split¬ 
ting between the two polarized emissions. This is a conse¬ 
quence of the sp-d interaction. At higher magnetic-fields this 
splitting is superimposed by opposed oscillations originated 
from the QW holes LLs oscillations. We observe that the non¬ 
linear behavior for the cr_ transition increases as the QW hole 
gas starts to be polarized reaching a maximum value when 
the hole gas is spin-polarized. As this polarization starts to 
decrease, the non-linear behavior for a- decreases. The op¬ 
posite behavior is observed for the cr+ transition. It shows a 
negative non-linear behavior which roughly follows the same 
dependence. As a consequence, if we look at the transition 
energy shift, (E ^~) shows a maximum oscillation at odd 
(even) filling factor. Combining all the effects, the non-linear 
splitting oscillates between 1.2 meV (B = 3.8 T) and a max¬ 
imum of 2 meV {B = 15 T). These values are in qualitative 
agreement and of the same order of the non-linear splitting 
observed experimentally. 

Our results suggest that the oscillations observed experi¬ 
mentally by Gazoto et al. [Ref. 14] are a consequence of a 
combined effect of the spin polarized hole gas, the Coulomb 
exchange interaction and the p-d exchange interaction. How¬ 
ever, in order to observe similar Mn dependence in the effect 
we consider the Mn ions closer to the QW as compared to 
their samples. 

A. Role of interactions and Mn position on the electronic 
structure 

We shall now investigate how the spin-dependent interac¬ 
tions and Mn position affects the electronic properties of the 
2DHG. We start our discussion considering the same sample 
parameters as discussed above, but turning off the sp-d inter¬ 
action, with Mn ions acting as non-magnetic acceptor impu¬ 
rity. This give us a clear picture of the effects induced by 
Coulomb exchange interaction on the hole gas, and conse¬ 
quently, by comparison with the previous section, we extract 
the effects of Mn on the hole gas. In the sequence, we change 
the Mn-doping position to 3 nm and turn-off the Coulomb ex¬ 
change interaction. This allows us to infer about the effects 
Mn position on the electronic structure. Here we will focus 
on the hole gas LLs and transition energy behavior with the 
magnetic field. The CB electronic states have the same quanti¬ 
tative behavior as in the case described in the previous section. 

Figure 4(a) illustrates the hh-LLs fan diagram in absence 



FIG. 4. (Color online)(a) Heavy-hole Landau Level fan diagram. 
Solid, dash-dot, dash-three-dots, dash-two-dots, and dashline repre¬ 
sents ^^, 4 , and Ep, respectively, (b) Transi¬ 

tion energy (right j/-axis) and non-linear energy shift (left j/-axis) as 
a function of the magnetic field. We considered La = 1 nm and 

Trhh(lh) _ T re _ ej 

Wd — ^sd - u. 


of the p-d exchange coupling. We first observe that the LLs 
always increase with the magnetic field independent of the 
hh-spin, differently of what was observed with the full Hamil¬ 
tonian for the same heterostructure, but in the presence of the 
hole-Mn spin interaction. The transition energies, shown in 
Fig. 4(b), for both polarizations, increases linearly with the 
magnetic field, as expected, and show oscillations that are di¬ 
rectly related with the LLs filling factor. A clear picture of the 
oscillations is given in the non-linear energy shift. For low- 
magnetic fields (B < 2T) the non-linear energy shift presents 
a week dependence with B. For higher-magnetic fields it os¬ 
cillates with LLs filling factor, with the maximum (minimum) 
shift at odd (even) filling factors. This result indicates that the 
dominant effect of the sp-d interactions at low magnetic fields 
with the consequent energy splitting between the different po¬ 
larizations while at higher magnetic fields the spin-polarized 
hole gas occupation is the dominant effect being responsible 
for the alternate oscillations observed. 

Following, we consider the Mn layer located at 3 nm from 
the QW/barrier interface. This is the nominal position of the 
Mn layer in the samples investigated in Ref. [14]. All the 
other parameters are the same as discussed before. We set 
now Vxc = 0 in the Hamiltonian, Eq. (2), to focus on the 
remaining sp-d interactions. Fig. 5(a) shows the holes LLs 
fan diagram. We observe the non-linear energy splitting be¬ 
tween the two polarizations at low magnetic fields. This split¬ 
ting, however, is considerable diminished as a consequence of 
the weaker overlap of the hole wave-functions with the Mn 
layer. We also observe an oscillatory behavior in the hh}^ 
LLs with the magnetic field which is associated to the QW LL 
filling factor. These oscillations, however, have a completely 
different qualitative and quantitative behavior in comparison 
with the results shown in the previously section. First of all. 
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we do not observe a qualitative dependence with the hole-spin 
in the LLs oscillations. Both fundamental states, hh^ ^ and 
hh^ oscillate following the same pattern with the magnetic 
field. Second, the value of these oscillations is significantly 
diminished. This can be better visualized in Fig. 5(b) where 
we plot the energy transition and the non-linear energy shift 
as a function of the magnetic field. The non-linear behavior is 
the same for both circularly polarized transitions, as expected 
from the hh LLs behavior. The tiny split observed in the total 
energy shift is of the order of 0.1 meV. These results con¬ 
firm the combined origin of the oscillations in the Coulomb 
exchange interaction and the p-d exchange interaction in the 
presence of the hole gas. They also demonstrate that the Mn 
ions have to be close to the QW states to observe a sizeable 
p-d coupling. 



FIG. 5. (Color online) (a) LL fan diagram for heavy-holes. The en¬ 
ergy states are, in order of increasing energy, labelled as ^ (solid), 
(dash-dot), (dash-three-dots), (dash two-dots) and 

i?F (dash-line) , respectively (b) Transition energy (right y-axis) and 
non-linear energy shift (left y-axis) as function of the magnetic field. 
We considered Ls = 3 nm and vxc = 0. 


B. Gate Voltage 

It is interesting to consider now the possibility to control 
the magnetic effects in the structure by the application of an 
electric field. This can be achieved through a gate voltage 
modifying the Fermi level in the structure and, therefore, the 
carrier distribution. Essentially, we will simple change the 
Fermi level at the surface by the expression Ep —>■ Eg/2 — Vg 
and will consider the value Vg = 0.71 eV which leads to an al¬ 
most flat band condition near the surface. We consider a value 
of 1 nm for the Mn layer spacing. All the other parameters are 
the same. 

Figure 6(a) shows the hh potential profile and the wave- 
functions of the occupied hh levels at (a) B=0 T and for (c) 
B=10 T, (b) the hh LLs fan diagram and (d) the hh-LLs charge 
concentration as a function of the magnetic field. We first ob¬ 


serve that there are three occupied hh states at B = 0 T. As 
the magnetic field increases, this picture changes completely. 
The Mn rich region becomes attractive for the spin up hhs 
and repulsive for spin down hhs, and the QW spin up levels 
are pushed towards this region. As a consequence, follow¬ 
ing the previous terminology, hh^ strongly overlaps with the 
Mn region. The same happens for hh^ although it is less pro¬ 
nounced. hh^ also presents a significant overlap with the Mn 
region, while /i/i| is mainly confined in the QW. The states 
hh^. and hh^ become almost fully confined in the C layer. 
All this evolution in the wave-functions produces a complex 
structure in the LLs fan diagram. 



FIG. 6. (Color online) (a) and (c) Self-consistent profile potential and 
heavy-hole wave-functions for B = 0 T and B = 10 T, respectively, 
(b) Heavy-hole LLs and (d) two-dimensional density of the LLs. The 
LLs and its density are, in order of increasing energy, labelled as 
(solid), (dash-dot), (dash-three-dots), ^ (dash 

two-dots) and Bp (dash-line), respectively. For XMn = 0.4 MLs, 
Vg = 0.71 eV and Ls = 1 nm. 


We focus now our analysis in the hh ground-state LLs, that 
is, /i/ig and hh^ At low magnetic fields, the two states 
show a strong non-linear behavior, with hh^ ^ showing a con¬ 
vex curvature while hh^ ^ shows a concave curvature with the 
magnetic field. This regime is entirely dominated by the p-d 
interaction. At higher magnetic fields, for B > 6 T, the oscil¬ 
lations with the LLs filling factors become more pronounced 
and dominate the features. However, in this case, we do not 
observe the difference in the magnetic field dependence from 
spin up and spin down with the odd and even filling factors 
which was originated from the hole gas exchange energy in 
the presence of the spin-polarization. Here, both states os¬ 
cillate in a similar way. The main reason for this behavior 
is in the way that the LLs associated to the QW are depop¬ 
ulated. This can be observed in Fig. 6(d). In particular, for 
H ~ 13 T, both LLs associated to spin up and down become 
depopulated almost at the same magnetic field. This prevents 
the spin-polarization of the QW hole gas. At the same time, 
the participation of the states hh'^ in the polarization of the 
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QW hole gas prevents a clear oscillation in this polarization. 
On the other way, the states /i/ij ^ have a strong influence of 
the Mn ions which dominates the magnetic held dependence. 
The oscillations are therefore mainly dominated by the charge 
transfer among the QW levels but not in the polarization of 
the hole gas. 
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FIG. 7. (Color online) (a) CB LL fan diagram . (b) Transition energy 
(right y-axis) and non-linear energy shift (left i/-axis) as function of 
the magnetic field. XMn = 0.4 MLs, Vg = 0.71 eV and Ls = 1 nm. 


Figure 7 shows (a) the CB LLs fan diagram and (b) the 
transition energies and the non-linear energy shift for this 
case. The CB LLs show the usual oscillation associated to 
the charge transfer between different hole gas reservoirs. The 
non-linear shift observes a rapid split at low magnetic fields 
due to the increasing of the sp-d interaction with the magnetic 
field. For B > 3T this splitting saturates and the oscillatory 
behavior dominates its features. As it was already mentioned 
for the hh LLs, we observe symmetric oscillations for the two 


circularly polarized transitions associated to the LLs filling 
factor. These results show that by applying a gate voltage we 
are able to qualitatively change the behaviour of the optical 
emission. 


IV. CONCLUDING REMARK S 

We presented here the results of calculations of the elec¬ 
tronic structure and optical emission energy for (In,Ga)As 
QWs with barriers doped with Mn and C acceptors with a fo¬ 
cus on the influence of the Mn impurities. Our results clearly 
show that for a significant magnetic effect on the optical emis¬ 
sion the Mn layer has to be near the QW interface, allowing a 
strong overlap between QW wave functions and Mn spin. We 
can separate two regimes for the optical emission effects. At 
low-magnetic fields, B < 5 T, the p-d exchange interaction 
dominates the effects while at higher magnetic fields the spin- 
polarization of the hole gas becomes the dominant effect. We 
also demonstrated that the non-linear effects can be controlled 
applying an external electric field. 

Our results explain the experimental observations by 
Gazoto et al. [Ref. 14]. The oscillatory magneto-shift ob¬ 
served in the magneto-photoluminescence can be interpreted 
as a combined effect of the exchange hole interaction in the 
presence of the spin-polarized hole gas enhanced by the p-d 
interaction. However, to observe a clear dependence with the 
Mn it is necessary to consider a Mn layer closer to the QW 
than the nominal value. 
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